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Abstract 



The production of a tt-pair and its decay into a 6 fermion final state of different flavours in e~^e~ 
anniliilation at centre of mass energies typical for linear colliders is analyzed in the framework 
of the Standard Model. The results of calculation based on exact matrix elements at the tree 
level and full 6 particle phase space are compared with a few different approximations. It is 
shown that the effects related to off-shellness of the tt-pair and background contributions are 
sizable both in the continuum and at the threshold. 



^Work supported in part by the Polish State Committee for Scientific Research (KBN) under contract No. 
2 P03B 004 18 and by European Commission's 5-th Framework contract HPRN-CT-2000-00149. 
^E-maih kolodzie@us.edu.pl 



1 Introduction 



The physical properties of the top quark directly measured at the Tevatron are in a very good 
agreement with those derived from the Standard Model (SM) analysis of the data collected 
at LEP and SLC However, as the top quark is the heaviest particle ever observed, with 
the mass close to the scale of electroweak symmetry breaking, the measurement of its Yukawa 
coupling may give hints towards better understanding of the electroweak symmetry breaking 
mechanism and observed fermion mass hierarchy. Should the effects of the physics beyond 
the SM be visible at the energy scale below 1 TeV, it is very likely that precise measurements 
of the top couplings to electroweak gauge bosons, or its electric and magnetic dipole moment 
show deviations from the corresponding SM values. The high-precision of measurements of the 
top quark properties and interactions can be by far the best reached at an e~^e~ collider which 
operates at a clean experimental environment. Therefore, such measurements are planned at 
TESLA [0] and will most certainly belong to the research program of any future e~^e~ collider 



It is clear, that in order to disentangle the possible new physical effects from physics of the 
SM, it is crucial to know the SM predictions for the top quark pair production and decay as 
precisely as possible. Due to the large mass and width, the top quark decays before toponium 
resonances can form and the predictions for 

e+e^ tt (1) 

can be obtained within the perturbative QCD. The predictions for reaction (|l|) in the threshold 
region were obtained in |Q and then improved by calculation of the next-to-next-to-leading 
order QCD corrections 0, and by including the effects of initial state radiation and beam- 
strahlung 0. The 0{aas) and 0{aa1) |§] corrections to the top decay into a W boson 
and a h quark are also known. In the continuum above the threshold, the QCD predictions 
for reaction (|l]) are known to order a1 and the electroweak (EW) corrections to one-loop 
order |T^, including the hard bremsstrahlung corrections ||TT|]. The QCD and EW corrections 



are large, typically of (9(10%). Order as [0 and a1 QCD, and EW corrections has been 
combined in 

As measurements of the top quark physical characteristics, in particular its static properties 
such as magnetic and electric dipole moments, will be performed at high energies, much above 
the tt threshold, it is crucial to know off-resonance background contributions to any specific 
6 fermion decay channel and to estimate the effects related to the off-shellness of the tt-pair. 
Therefore, in the present note, instead of considering production of the top quark pair (|l|) and 
its subsequent decay into a specific 6 fermion final state, the 6 fermion reactions of the form 

e+e- ^ bfJibhh, (2) 

where fi = i'^,i't, u, c, /2 = fi^ ,d, s, and f[, f!^ are the corresponding weak isospin partners 
of /i, /2, fi = fi~ ,T~ ,d, s, f2 = i^^i,i^T, u, c, are studied in the lowest order of SM. For the sake 
of simplicity, it is assumed that the actual values of /i and are different from each other, 
and that neither f[ nor /2 is an electron. The results for reaction (Q) are compared with the 
results obtained in a few different approximations: the double resonance approximation for 
W bosons 

e+e- ^ bW+*bW-* ^ bfJibf2h (3) 



where only those 61 Feynman diagrams are taken into account which contribute to e+e^ — 
bW~^bW~ and the W bosons are considered as being off-mass-shell, the double resonance 
approximation for a t- and t-quark 

e+e- ^ ft* ^ bfJibhh (4) 

with only two 'signal' diagrams contributing and, finally, with 3 different narrow width ap- 
proximations: for the W bosons, top and antitop quarks, and a single top quark |l5l . 



A similar analysis of the 6 fermion processes relevant for a ti production in e~^e~ annihilation 
have been performed in refs. ||T^, where semileptonic channels of reaction have been 
studied, and in ref . JlSl , where purely hadronic channels of @ have been analysed. Moreover, 
irreducible QCD background to top searches in semileptonic channels of @ has been discussed 
in [|19]. The novelty of the present work, besides the more detailed discussion of the different 
approximations listed above, consists in taking into account both the electroweak and QCD 
lowest order contributions. Moreover, as light fermion masses are not neglected, the cross 
sections are calculated without any kinematical cuts. 

Basics of the calculation are described in the next section. Numerical results are presented 
and discussed in Section 3 and, finally, in Section 4, the concluding remarks are given. 



2 Calculation 



The calculation of matrix elements of reaction (0) is based on the complete set of the Feynman 
diagrams at the tree level of SM. The number of diagrams which contribute to (0) in the unitary 
gauge, neglecting the Higgs boson coupling to fermions hghter than a b quark, amounts to 
201 for semi leptonic final states, which contain two different charged leptons, and to 333 for 
purely hadronic final states, with different quark flavours. The necessary matrix elements are 
calculated with the method proposed in and further developed in |^ . As in , fermion 



masses are kept nonzero in the matrix elements and in the kinematics. The constant widths of 
unstable particles, the massive electroweak vector bosons, the Higgs boson and the top quark, 
are introduced through the complex mass parameters: 

My = my- imvTv, V = W,Z, = rriH - imnTH, Mt = rrit - iTt/2, (5) 

which replace masses in the corresponding propagators, both in the s- and t-channel Feynman 
diagrams. 



Propagators of a photon and a gluon are taken in the Feynman gauge. 



The 6 particle phase space of reaction _ 

^4 L_ ^ - rr 

.=3 /.=3(2^)'2E. 



where the energies and momenta of the initial state particles of reaction has been numbered 
from 1 to 2, and those of the finale state particles from 3 to 8, is parametrized in three different 



ways 

]14 r 1 //o„\14 



d Lips = l/(27r) dPS2 (s, S345, Says) dPS2 {ss45,m^, S45) dPS2 [sers^mQ, s^g 



X 



dPS2 (S45, m^, dP^z (syg, dS345dS678dS45dS78, (8) 



d^^Lips = l/(27r)^^dPS'2 (s, S34, ssers) dPS'2 (sgers, S56, 573) dPS'2 (534, mj, mfj 

X dPS'2 (S56, "^e) dP5'2 (S78, m^, mg) dS34dS5678dS56dS78. (9) 

and 

d^^LipS = 1/ (27r)^^dPS'2 (s, mg, 545673) dPS'2 (S45678, S45, ^678) dPS'2 (S678, "^6, S78) 



X 



dPS2 (s45, ^4, mfj dPS2 (S78, my, mg) dS45678dS45dS678dS78. (10) 

In Eqs. (|8HT0|), Sijfc... = (Pi+Pj +Pk + ■■■Y, hj, k = 3, 8, and dPS'2 (s, s', s") is a two particle 
(subsystem) phase space element defined by 



dP^2 (s, s', s") =6^{p-p'- p") ^^^V _ ^dfi', (11) 
^ ' ' ^ ^ ^ ^ ' 2E' 2E" A^/^ ' ^ ^ 

where p ' is the momentum and Q' is the solid angle of one of the particles (subsystems) in 
the relative centre of mass system, p' + p" = 0. Using the rotational symmetry with respect 
to the beam line, an integration over one azimuthal angle in the c.m.s. becomes trivial. This 
reduces the number of necessary integrations to be performed. 

Parametrization (P) is most suitable for integrating the dominant tt resonance contributions 
of Eq. (§). Parametrization (P) covers best contributions corresponding to the double W 
resonance approximation of Eq. (|^), whereas parametrization (p!OD covers other 'background' 
contributions to reaction (^. Parametrizations (|8|-[lOD are used with different permutations 



of external particle momenta and with different mappings which take into account the Breit- 
Wigner shape of the W, Z, Higgs and top quark resonances as well as the exchange of a 
massless photon or gluon. For a given final state and c.m.s. energy, altogether about 60 
kinematical channels are sampled in order to find the dominant channels which contribute 
more than 0.1% to the total cross section. Those dominant kinematical channels are then 
used in a multichannel Monte Carlo (MC) integration routine. 

The phase space integration is simplified in the narrow width approximations. The cross 
section of reaction (|]) in the narrow W width approximation is given by 

'TbwHw- = ^(e+e- ^ bWnW-) T {w+ ^ F {W" ^ /2/^) /F^^. (12) 

Similarly, the cross section of reaction (^) in the narrow width approximation for the top and 
antitop reads 

au = a(e+e- ^ tt) T [t ^ hfj[) F (t ^ 6/2/^) /F,^. (13) 

Finally, in the approximation where only the top quark is put on its mass shell, the cross 
section is given by 

<ytlM!, = ^(e+e- ^ tbf2h) T{t~^ hhf[) /F^. (14) 

There are 7 integrations which have to be performed numerically in order to obtain total cross 
sections in approximations (|12|) and ( p!^ and only one integration in the case of approximation 



T3|) . All numerical integrations in the present work are performed with VEGAS p2[ . 



3 Numerical results 



In this section, numerical results for the total and a few differential cross sections of reac- 
tion (0) are presented. They are compared with the corresponding results obtained within 
approximations (D, (|) and (|12|-[l|). 



The SM electroweak physical parameters are defined in terms of the gauge boson masses and 
widths, the top mass and the Fermi coupling constant. The actual values of the parameters 
are taken from 

mw = 80.419 GeV, Tw = 2.12 GeV, mz = 91.1882 GeV, Tz = 2.4952 GeV, 

mt = 174.3 GeV, = 1.16639 x 10"^ GeV'^. (15) 

The Higgs boson mass is assumed to be rnn = 115 GeV and the Higgs width is calculated 
according to the lowest order of SM resulting in Th = 4.9657 MeV. The top quark width is 
taken to be = 1.5 GeV. 

The SM electroweak coupling constants are given in terms of the electric charge ew = 
(Anaw)^^'^ and electroweak mixing parameter sin^6'iy with 

aw = V^G^m'^ sin^ 6^/ tt, sin^ 6iv = I — m^/ m|, (16) 



where m^y and mz are physical masses of the and boson specified in Eq. ([I5|). This 
kind of parametrization, together with substitutions of Eq. (^, is usually referred to as the 
'fixed width scheme' (FWS). The strong coupling constant is given by Qs = {Anas^Mz))^^"^ , 
with as{Mz) = 0.1185. 

It is also possible to perform computations with the complex electroweak mixing parameter 

sin2^^ = l-M^/M|, (17) 

with and Mf defined in Eq. (^. This kind of parametrization is called the 'complex- 
mass scheme' (CMS) CMS has the advantage that it preserves the SU{2) x f/(l) Ward 
identities [p3[] . 

For the sake of definiteness, other fermion masses used in the calculation are listed below [|l| : 
me = 0.510998902 MeV, m^ = 105.658357 MeV, m^ = 1777.03 MeV, 
mu = 5 MeV, = 9 MeV, m, = 150 MeV, m^ = 1.3 GeV, mb = 4.4 GeV. (18) 
The Cabibo-Kobayashi-Maskawa mixing is neglected. 

Matrix elements of e^e^ hW^hW~ and e^e~ —>■ ti have been checked against MADGRAPH 
p5| showing an agreement up to 13-16 decimals. As the version of MADGRAPH used in the 



comparisons is not applicable to processes with 6 particles in the final state, it has been 
not possible to compare directly matrix elements of reactions (|^) with those generated by 
MADGRAPH. Instead of that, matrix elements of different 'subprocesses' of (§), namely e+e^ 
bfif[bW~, e~^e~ bW^bf2f2 and e~^e~ fififif2^, have been compared successfully. The 
multichannel phase space generation routine has been checked by comparing normalization of 



different cliannels against eacli otlier and testing energy-momentum conservation and on-mass- 
sliell relations. For several total cross sections, ttie numerical integration lias been performed 
with different parametrizations of the phase space and the results have been stable within one 
standard deviation. 

The standard deviation of the multichannel integration routine is obtained as a sum of the 
standard deviations calculated by VEGAS for individual channels. This gives a more conserva- 
tive estimate of the integration error than for example adding up partial errors in quadrature. 

Another test is a comparison with existing calculations. Results for total cross sections of 
e'^e" bUfj^fi'^bdu in the lowest order SM are compared with the results of jl^ in Table 1. As in 



ITj , only the pure electroweak diagrams and the two tt signal diagrams are taken into account 



in (JaiiEW and at*t*, respectively. For the sake of comparison, the physical parameters of |]T7 



have been used, i.e. mz = 91.187 GeV, Tz = 2.49 GeV, mw = 80.22 GeV, Tz = 2.052 GeV, 
mt = 174 GeV, Tt = 1.558 GeV, = 4.1 GeV, m„ = 2 MeV and = 5 MeV. The 
electroweak mixing parameter is defined as in Eq. (|16D and aw = 1/128.07 is used at the same 
time. As the values of the Higgs boson mass and width used in the calculation are not quoted 
in [O], ttlh = 115 GeV and the lowest order SM value Th = 4.3977 MeV, corresponding to 



parameters of [|T3 , have been used in Table 1 . Another source of ambiguity in the comparison 



is the treatment of the finite widths of unstable particles which is not explicitly described in 



171 . Therefore the prescription of Eqs. (||) and have been adopted. The results for ciaii ew 
and at*t* are shown in columns 2 and 3, whereas the corresponding results of the present 
work are shown in columns 4 and 5. The results for at*t* agree nicely within the uncertainties 
quoted in parenthesis. The agreement is still nice for the complete electroweak cross sections 
CaiiEW above the ti threshold. Below the threshold, at = 340 GeV, there is a substantial 



relative discrepancy between o"aiiEW of [1^ and that of the present work. It is amazing that 



the results for all the approximated cross sections listed in Table 2 of [0 agree with the 
present work, also at ^/s = 340 GeV. It is difficult to state definitely what the actual reason 
for this discrepancy is. However, most probably it is the Higgs boson contribution, and in 
particular the Higgs-strahlung 'subrocess' e^e~ ZH with the Higgs boson decaying into 
a virtual Vr"'"iy~-pair that is responsible for it. The results of the present work without the 
Higgs contribution are shown in the last column of Table 1. They nicely agree with the results 
of which contain the Higgs with its mass and width not being specified. 



Unfortunately, a similar detailed comparison with results of []18| is not possible, as the authors 
do not specify numerical values of the physical parameters used in their computations. As the 
cross section of reaction (0) at tree level is of 0{a^,), it is very sensitive to the choice of initial 
parameters. Although it is meaningless to perform any quantitative comparison, the results 



of the present work are in a qualitative agreement with those of ref. |T^ which will be shown 
later. A detailed quantitative comparison with Accomando, Ballestrero and Pizzio |T^ is also 
not possible as the authors include some radiative effects in most of their results. A meaningful 
comparison could in principle be performed for the Born cross sections of e~^e~ — > hv^^'^hdu 
corresponding to the ti signal and background at y/s = 500 GeV. With cuts of fl^ and the 
physical parameters of the present work, one obtains 17.895(9) fb and 1.25(2) fb for the signal 
and background, respectively. The result for the signal cross section differs from that of [|1^ by 
about 2% while the relative difference between the background cross sections is much bigger, 
probably because there is no gluon exchange contribution included in the Born background 



cross section of |jT6| 



Lowest order SM total cross sections of the semileptonic channel e~^e —>■ hv^^^hdu of reaction 



Table 1: Comparison of the lowest order SM total cross sections of e^e~ bv^fi^bdu of |17] and 
present work. Results of |jl^ obtained with a complete set of the electroweak diagrams, (XaiiEWi and 
two tt signal diagrams, cT^*f*, are shown in columns 2 and 3, whereas the corresponding results of the 



present work are shown in columns 4 and 5. Here the parameters of [17| are used. All cross sections 
are in fb. The number in parenthesis show the uncertainty of the last decimal. 



(GeV) 


F. Yuasa, et. al. [17 


Present Work 


0"all EW 




0"all EW 




no Higgs 
^all EW 


340 
350 
360 
380 
500 


0.687(2) 
6.45(1) 
14.97(2) 
21.42(4) 
22.32(4) 


0.4462(3) 
6.187(4) 
14.63(1) 
21.00(1) 
21.30(1) 


0769(10) 
6.59(2) 
15.03(5) 
21.48(8) 
22.55(4) 


0.4455(4) 
6.175(4) 
14.623(9) 
20.99(1) 
21.27(1) 


0.689(2) 

6.45(1) 
14.97(3) 
21.49(5) 
22.32(5) 



(0) at different c.m.s. energies typical for future linear colliders are shown in Table 2. The 
complete lowest order result a, the approximation of Eq. (^ 0"feTy+*6VK-*! the narrow W width 
approximation of Eq. ( |12|) cr^vK+feVK-! the approximation of Eq. (|) a ft*, the narrow width 
approximation of Eq. (p!3D for a top and an antitop quark and the narrow width approx- 
imation for a top quark of Eq. (0) cr^g^^ have been all obtained in FWS. The SM tree level 
analytic expression for the partial widths of the W boson and the experimental value of total 
W width Tw have been used in Eq. (p!2D. Similarly, the SM tree level analytic expression for 
the partial widths of the t quark in the zero fermion mass approximation and the total 



top width Tt = 1.5GeV have been used in Eqs. (|^) and ([Tj). The use of these values of Tw 
and Tt in Eqs. (|T2|-|Iip is preferred in the comparison because the same values have been used 
in substitutions of Eq. (|^). In Table 2, the numbers in parenthesis are standard deviations of 
the MC integration, which show an uncertainty of the last decimal. 

A cross section of the 6 fermion reaction e+e~ bu^^^bdu is nonzero already below the tt-pair 
production threshold. It is the single top (antitop) resonance and nonresonant background 
contributions which are responsible for that effect. Whether this background may affect phys- 
ical observables in the threshold region, such as the top invariant mass distribution or angular 
distributions of the final state quarks or leptons, will be discussed later. Close to threshold, 
at ^/s = 360 GeV, the relative difference between a and the narrow width approximation 
att is about -1.5%, whereas in the continuum the difference becomes bigger, as relevant as 
radiative corrections, amounting to 7% at a/s = 500 GeV and 19% at y/s = 800 GeV. At 
higher energies, the difference between a and att becomes so large that approximation (|13]) 
does not make sense any more. Comparison of approximated results ai,w+*bw-* ^md ai,w+bw- 
with the complete result a shows that approximations of Eqs. (|^) and (|I2]) are relatively much 
better in a wide range of the c.m.s. energy from 360 GeV to 2 TeV. 

The pure off-shellness effects of the tt-pair can be regarded as the difference between approx- 
imations at*t* of Eq. (H) and att of Eq. ([T3|) . They are plotted in Fig. 1 as a function of the 
c.m.s. energy. The two plots in Fig. 1 show a similar behaviour with c.m.s. energy as those in 
Fig. 4 of F. Gangemi et. al. |T^. A naive multiplication of the results plotted in Fig. 1 and 



Table 2: Lowest order SM total cross sections of e'^e" bi'f^fi^bdu in fb at different c.m.s. energies: 
the complete lowest order result a, the approximation of Eq. (|3|) aiji^+*i^r-* , the narrow W width 
approximation of Eq. ( [T^ ) a^^^r+iw-i the approximation of Eq. (Q) a^.j. , the narrow width approxi- 
mation of Eq. ([T3|) for a top and an antitop quark cr^f and the narrow width approximation for top 
quark of Eq. (|T^ ) cr^i^fi. The number in parenthesis show the uncertainty of the last decimal. 



(GeV) 


a 


'^bW+*bW-* 


'^bW+bW- 




(ytt 


'^tbdu 


340 


1.162(7) 


0.681(6) 


0.671(1) 


0.3521(2) 




0.2546(3) 


360 


13.64(2) 


13.224(8) 


13.618(8) 


12.79(1) 


13.875 


13.42(1) 


500 


20.48(9) 


20.17(1) 


20.79(1) 


19.06(1) 


19.223 


19.51(3) 


800 


10.61(4) 


10.46(3) 


10.75(1) 


9.181(5) 


8.918 


9.47(1) 


1000 


7.35(4) 


7.33(4) 


7.54(1) 


6.171(4) 


5.862 


6.390(7) 


2000 


2.43(2) 


2.48(3) 


2.48(1) 


1.847(2) 


1.510 


1.822(2) 



the results for at*t* of Table 2 by a factor 12, corresponding to the different colour factor and 
the sum over 4 different hadronic channels, gives nice agreement with the signal cross section 



plotted in Fig. 4 of 1 18 
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Figure 1. Total cross sections of e~^e~ bu^jjL'^bdu in approximations of Eqs. 

(P^) (dashed line) as functions of the c.m.s. energy. 



(solid line) and 



Lowest order SM total cross sections for different channels of (^ at c.m.s. energies typical 
for TESLA are compared in Table 3. The cross sections of e"'"e^ hcsbdu is about 3 times 



Table 3: Lowest order SM total cross sections in fb for different top production channels at c.m.s. 
energies typical for TESLA. The number in parenthesis show the uncertainty of the last decimal. 



(GeV) 






e+e bcsbdu 


360 


4.36(1) 


13.65(4) 


42.1(2) 


500 


6.70(2) 


20.48(9) 


62.2(2) 


800 


3.43(2) 


10.61(4) 


32.1(1) 



bigger than the cross section of e^e^ hv^jjL^hdu, which in turn is about 3 times bigger than 
the cross section of e+e~ hv^^^hr^Vr- This reflects the relative numbers of colour degrees 
of freedom. Small deviations of the relative factors from 3 result from the gluon exchange 
contributions, which are absent for e^e~ — >■ hv^^'^hr^Vr and are different for e^e~ hv^^'^hdu 
and e+e~ bcsbdu. The errors given in parenthesis have been obtained in the same way and 
have the same meaning as those of Tables 1 and 2. 
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^du (GeV) 



Figure 2. Differential cross sections da/dmf of e~^e~ bu^fi~^bdu at ^/s = 360 GeV versus the 
invariant mass of the t quark reconstructed from the bdu system (solid histogram) and from the 
bW~ system in the narrow W width approximation (dotted histogram). 



How the background nonresonant contributions affect differential cross sections of (^ is illus- 
trated in Figs. 2 and 3. In Fig. 2, the differential cross sections da/dm^^- of e~^e~ bu^fi'^bdu 
at a/s = 360 GeV are plotted versus the invariant mass of the i quark reconstructed from the 
bdu system. The three histograms: solid, corresponding to the complete lowest order result. 



dotted, corresponding to the approximation of Eq. (|^), and dashed one, representing the ti 
signal (^, shown in Fig. 2, are almost indistinguishable. This means that approximation (|[) is 
satisfactory and the background contributions coming from the single top (antitop) resonance 
and the nonresonant Feynman diagrams is negligible in this case. 

The differential cross sections dcr/dcos6' of e+e^ hu^^^hdu at y/s = 360 GeV are plotted 
versus cosine of the (up going curves) and d (down going curves) angle with respect to 
the positron beam in Fig. 3. The angular distributions obtained with the complete set of tree 
level Feynman diagrams differ substantially from the distributions based on approximations of 
Eqs. (1^) and (^. The final state muon fi~^ (down quark d) goes more preferably in the direction 
of initial positron (electron) than it would follow from the approximated distributions based 
on Eqs. (|) and (|). 
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Figure 3. Differential cross sections of e'^e — > bv^fi'^bdu at ^/s = 360 GeV versus cosine of the 
(up going curve) and d (down going curve) angle with respect to the positron beam. 

4 Summary and Outlook 

The production of a tt-pair and its decay into a 6 fermion final state of different flavours in 
e~^e~ annihilation at energies typical for linear colliders has been analyzed in the framework 
of the SM. The results of calculation based on exact matrix elements at the tree level and full 
6 particle phase space have been compared with the results obtained within a few different 
approximations: the double resonant approximations for the W bosons (^ and for the top 
and antitop quarks (^, the narrow width approximation for the W bosons (|12|), the narrow 
width approximation for the t- and t-quark (|l^) and for the t-quark quark only ([I^). 



It has been shown that the effects related to the off-shellness of the tt-pair and to presence of 
background contributions to cross sections of six fermion reactions (0) are quite substantial 
They are at the level of a few per cent already in the ti threshold region. In the continuum, 
at higher energies, the effects become quite sizable, reaching about 20% at ^/s = 800 GeV. 
Therefore, for achieving the desired precision level in the analysis of experimental data from 
linear colliders, it is mandatory to include them in theoretical predictions together with ra- 
diative corrections. The inclusion of the latter should reduce the dependence on the choice 
of initial parameters mentioned in Section 3 in the context of comparisons with the existing 
calculations. 
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